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Chapter  1 


Introduction 

In  the  past  the  lack  of  conventional  meteorological  observations,  especially  in  the 
less  accessible  regions  of  the  world  the  oceans  and  sparsely  populated  regions), 
hindered  a  lot  of  the  study  of  weather  systems  in  the  mentioned  areas.  Routine  analyses 
and  forecasts  had  to  rely  on  small  amounts  of  observational  data. 

The  analysis  and  forecast  of  meteorological  (especially  mesoscale)  phenomena  was 
also  hindered  at  intervals  between  the  synoptic  observation  times  due  to  the  non-availability 
of  conventional  data.  Conventional  data  means  the  data  of  ground-based  meteorological 
observation  instniments  such  as  thermometers,  and  upper-air  balloons,  or  radiosondes. 
Nevertheless  the  availability  of  observations  in  between  main  synoptic  observation  times  is 
critical  for  the  forecasting  of  meteorological  hazards  in  many  Helds,  such  as  aviation. 
Knowledge  early  in  time  of  hazardous  events  like  squall  lines  or  severe  storms  can  provide 
significant  warning  potential  and  thus  may  save  lives. 

With  the  availability  of  polar  orbiting  siUellites  like  the  NOAA/TIROS-N  (Tele¬ 
vision  and  Infrared  Operational  Satellite)  the  quantitative  analysis  of  the  three-dimensional 
state  of  the  atmosphere  is  possible  at  times  between  synoptic  observations  over  all  regions 
of  the  globe.  Due  to  the  orbit  characteristics  (sun-synchronous,  orbit  period  102  min.)  the 
whole  globe  is  covered  two  times  a  day  by  each  satellite.  With  a  swath  width  of  about 
2500  km  there  is  a  more  frequent  area  overlap  at  higher  latitudes. 

It  is  possible  to  gain  information  on  mesoscale  phenomena  at  the  resolution  of  the 
satellite  instruments.  The  TIROS  Operational  Vertical  Sounder  (TOYS)  instrument  pro¬ 
vides  data  that  make  it  possible  to  determine  the  three-dimensional  structure  of  the 
atmosphere  at  a  scale  of  about  100  km  and  largo*  in  horizontal  space  as  well  as  1-3  hours  in 
time,  which  is  the  definition  of  the  Meso-^  scale  (after  Orlanski,  1975).  Quantitative 
information  on  atmospheric  temperature  fields  can  be  obtained  with  considerable  accuracy 
(LeMarshall,  1988).  Recently,  there  have  been  efforts  to  estimate  the  atmospheric  humidity 
with  TOVS-derived  data  (Husson  et  al.,  1991).  Up  to  now  there  are  not  many  examples 
available  on  the  study  of  mesoscale  water  v^or  structure  in  the  atmosphere  using  only 
satellite  information. 

Several  studies  of  mesoscale  phenomena  have  been  carried  out  in  the  past,  covering 
the  Antarctic  (Heinemann,  1990)  the  arctic  region,  and  polar  lows  (Claud  et  al.,  1990, 
Prangsma,  1989)  as  well  as  mid-latitude  cloud  clusters  (Klaes,  1991). 
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The  purpose  of  this  research  is  to  demonstrate  the  capabilities  of  two  TOYS 
retrieval  algorithms.  The  algorithms  used  here  are  the  Improved  Initialization  Inversion 
(31)  (Chedin  et  al.,  1985)  as  well  as  the  International  TOYS  Retrieval  Package  (TIPP) 
methods  (Smith  et  al.,  1985).  In  this  report  the  application  of  the  TOYS  retrieval 
algorithms  “SF’  and  ITPP  is  demonstrated  with  the  example  of  a  thunderstorm  situation  in 
summer  over  central  Europe.  It  occurred  on  12-13  August  1990.  One  significant  event  is 
the  development  of  a  large  mesoscale  cloud  cluster  over  eastern  France  and  Southwest 
Germany. 

In  Section  2,  there  is  an  overview  of  the  TOYS  instruments,  including  the  HIRS/2 
High  Resolution  Infrared  Sounder  and  the  Microwave  Sounding  Unit  (MSU),  both  of 
which  are  used  by  the  retrieval  algorithms.  The  two  retrieval  algorithms  will  be  discussed 
in  Chapter  3.  Chapter  4  describes  the  synoptic  situation  over  central  Europe  on  12-13 
August  Chapter  5  shows  some  retrieval  results  obtained  by  the  two  algorithms  for  the 
case  study.  Finally,  Chapter  6  gives  a  sununary  of  the  possible  comparisons  with  other 
data  sources  including  the  conventional  meteorological  data. 


2 


Chapter  2 


The  TOYS  Sensors 

The  TOVS-instrument  aboard  the  NOAA  satellites  is  comprised  of  three  passive 
sounding  instruments  (Smith  et  aL,  1979).  The  HIRS/2  is  a  radiometer  with  20  channels, 
the  MSU  is  a  passive  microwave  radiometer  with  four  channels,  finally  the  Stratospheric 
Sounding  Unit  (SSU)  is  an  infrared  radiometer  with  three  channels,  based  on  the  principle 
of  pressure  modulation.  The  TOYS  instrument  characteristics  are  contained  in  Planet 
(1988).  Therein  is  given  detailed  information  on  the  calibration  characteristics  and  the  data 
formats.  The  following  subsections  list  the  radiative  sensing  characteristics  of  each  unit 

2.1  High  Resolution  Infrared  Radiation  Sounder, 

Model  2  (HIRS/2) 

The  HIRS/2  is  a  further  developed  version  of  the  HIRS/1  which  flew  on  the 
NIMBUS-6  satellite.  It  measures  upwelling  radiation  in  19  regions  of  the  infrared  spec¬ 
trum  and  one  region  of  the  visible  spectrum.  The  ground  resolution,  i.e.  the  foo^rint  of 
the  HlRS/2  is  17.4  km  in  diameter  at  nadir.  There  are  56  scanning  steps  in  one  scan-line, 
one  scan  takes  100  milliseconds,  and  the  complete  scan-line  is  scanned  in  6.4  seconds. 
The  scan  covers  an  angle  of  ±  49.S  deg,  or  ±  1125  km  left  and  right  from  the  subsatellite 
point  for  a  swath-width  of  2250  km.  Due  to  the  curvature  of  the  earth  the  footprint  at  the 
edge  of  scan  is  different  from  the  one  at  nadir:  58.5  km  cross-track  by  29.9  km  along 
track.  That  means  that  there  is  an  elliptical  pattern  of  the  HlRS/2  spots  off  the  subsatellite 
track.  Figure  1  shows  the  ground  pattern  of  the  HIRS/2  instrument 

After  every  40  HIRS/2  scan  lines  there  is  an  in-flight  calibration  of  the  instrument 
during  which  there  are  no  retrieval  data  available.  This  means  a  g^  of  tiiree  lines  for  each 
calibration  cycle. 

Table  1  shows  the  characteristics  of  the  HIRS/2  channels.  There  are  seven  chan¬ 
nels  (1-7)  in  the  15  pm  region,  mainly  used  for  temperature  sounding.  Channel  8  at  1 1.10 
pm  is  a  window  channel  and  serves  for  cloud  detection  and  surface  temperature  estimation. 
Channel  9  at  9.7  pm  is  the  ozone  channel.  Total  ozone  concentration  may  be  obtained  by 
using  its  data.  Channels  10-12  are  used  for  water  vapor  retrieval.  They  are  positioned 
from  8.3  pm  to  6.7  pm.  Further  information  on  the  temperature  profile  of  the  atmosphere 
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is  retrieved  by  using  the  information  of  HIRS/2-channels  13'17,  which  are  positioned 
around  4.3  pm.  They  provide  a  better  sensitivity  to  wanner  regions  of  the  atmosphere  than 
do  the  13  iim  channels.  They  are  also  less  sensitive  to  clouds  than  the  IS  ^m  channels. 
There  are  two  more  window  channels  with  central  wavelengths  at  4.0  and  3.7  ^im.  They 
too  can  be  used  to  estimate  the  surface  temperature  and  to  detect  cloud  presence.  Finally, 
there  is  a  channel  in  the  visible  part  of  the  spectrum.  Channel  20  at  0.7  p.m  is  used  during 
daytime  for  cloud  detection  and  thus  to  help  define  clear  fields-of  view. 

Figure  2  .  .ows  the  weighting  functions  for  the  HIRS/2  channels.  Each  peak 
represents  the  region  of  the  atmosphere  from  which  a  channel  receives  its  major  energy 
contribution. 

2.2  Microwave  Sounding  Unit  (MSU) 

The  MSU  is  a  passive  microwave  Dicke  radiometer  (see  e.g.  Vowinkel,  1988).  Its 
four  channels  are  positioned  in  the  50  GHz  region  of  the  electromagnetic  spectrum  (the 
5.5  mm  oxygen  band). 

The  horizontal  scan  pattern  is  109.3  km  at  nadir.  There  are  1 1  scans  at  a  scan  time 
of  25.6  seconds  per  scan.  As  shown  in  Figure  1,  the  cross-track  scan  angle  is  ±  47.35 
deg,  which  means  a  swath  of  2  x  1 173  km,  or  2346  km.  The  field-of-view  pattern  at  the 
end  of  a  scan  is  323.1  km  cross-track  and  178.8  km  along  track. 

Channels  2-4  are  used  for  temperature  sounding  through  clouds.  Channel  1  at 
50.31  GHz  is  a  window  channel  and  is  used  to  determine  cloud  attenuation  and  surface 
emissivity  parameters.  Table  2  displays  the  characteristics  of  the  MSU  channels.  Figure  3 
shows  the  weighting  functions  of  the  MSU  channels. 

2.3  Stratospheric  Sounding  Unit  (SSU) 

The  SSU  measures  the  absorption  of  radiation  by  carbon  dioxide  contained  within 
its  sensors.  The  pressure  of  the  carbon  dioxide  in  the  gas  cell  determines  the  altitude  of  the 
weighting  function  peaks  in  the  atmosphere.  Table  3  gives  the  characteristics  of  the  SSU 
instrument  channels.  The  SSU  will  not  be  presented  further,  because  it  is  not  used  for  the 
retrievals  discussed  in  this  study. 


4 


Chapter  3 


Retrieval  Algorithms 

It  is  well  known  that  the  three-dimensional  structure  of  the  atmosphere  may  be 
deduced  in  terms  of  temperature  and  moisture  from  measurements  of  the  earth’s  and 
clouds’  radiance,  measured  by  satellite  or  aircraft-based  multichannel  passive  sounding 
instruments.  This  is  done  by  a  so-called  “inversion”  procedure,  which  in  the  present  case 
means  die  inversion  of  die  radiative  transfer  equation. 

There  are  several  problems  involved  in  die  inversion  procedure.  One  is  the  under- 
determination  of  the  problem,  making  it  necessary  to  obtain  additional  information  from 
some  other  source.  Another  is  the  non-uniqueness  of  the  inversion  of  the  radiative  transfer 
equation  -  that  means  given  a  state  of  the  atmosphere  in  terms  of  temparature  and  moisture 
profiles,  and  information  on  the  important  gaseous  constituents  of  the  atmosphere,  then  a 
unique  solution  for  the  vertical  distribution  of  radiances  may  be  obtained.  On  the  other 
hand,  inverting  the  radiative  transfer  equation  by  having  as  input  a  vertical  profile  of 
radiances  yields  a  multiplicity  of  solutions  for  the  three-dimensional  state  of  the  atmosphere 
in  terms  of  temperature  and  moisture,  depending  on  the  respective  situation.  To  solve  the 
problem,  a  so-called  *Tirst  guess”  is  necessary  to  start  the  solution  process.  The  solution  is 
very  strongly  dependent  on  the  choice  of  the  initial  guess,  especially  in  regions  (like  in 
cloud  areas)  where  less  information  is  obtainable  (Susskind  and  Rosenberg,  1980).  A 
detailed  discussion  of  the  inversion  problem  is  given  by  Twomey  (1977). 

Different  software  packages  are  available  to  derive  atmospheric  parameters  from 
TOVs  data.  These  packages  make  use  in  varying  degrees  of  two  basic  types  of  retrieval 
methods:  1)  Statistical  retrievals  begin  with  the  measured  radiances  and  apply  regression 
schemes  on  them  to  obtain  the  desired  atmosphoic  parameters  such  as  temperature  or  water 
vapor,  statistical  retrievals  are  fast,  making  them  for  a  long  time  the  only  candidates  for 
operational  application.  The  disadvantages  are  the  large  amount  of  reference  situations 
needed  to  estimate  the  regression  coefficients,  and  following  from  this  there  are  misleading 
results  in  atmospheric  situations  deviating  from  standard  situations.  2)  Physical  retrievals 
solve  the  radiative  transfer  equation  using  an  iteration  process  until  a  sufficient  good 
solution  matching  the  actual  observations  is  reached.  Physical  retrievals  require  the 
computation  of  transmittances  and  radiances  at  each  iteration.  They  also  need  a  “first 
guess”  which  is  the  more  important  factor,  since  when  it  is  badly  chosen  it  influences  the 
result  in  a  negative  sense. 
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A  physical  scheme  is  presently  used  operationally  by  NOAA/NESDIS  (Heming  et 
al.,  1986).  It  has  replaced  a  statistical  scheme  (Reale  et  al.,  1986).  One  other  widely  used 
scheme  is  the  “International  TOYS  Processing  (ITPP)  Package,”  developed  and  maintained 
at  the  University  of  Wisconsin,  Madison  (Smith  and  Woolf,  1976;  Smith  et  al.,  1979; 
Smith  et  al.,  1985).  ITPP  is  a  physical  schen^  that  allows  the  simultaneous  retrieval  of 
surface  temperature,  the  vertical  temperature  profile,  and  the  vertical  moisture  profile. 

The  Improved  Initialization  Inversion  (“ST*)  algorithm  (Chedin  and  Scott,  1984; 
Chedin  et  al.,  1985;  Chedin  and  Scott,  1985),  developed  and  maintained  at  Laboratoire  de 
Meteorologie  Dynamique  (LMD)  at  Palaiseau,  France,  is  a  physico-statistical  method.  It 
relies  upon  a  priori  knowledge  of  the  atmospheric  structure  and  a  brightness  temperature 
(Bl^  for  each  channel  used  of  the  TOVS-instrument 

One  other  method  is  the  differential  retrieval,  a  direct  method  which  has  not  been 
^lied  to  operational  data  yet  (King,  1990;  Ou  et  al.,  1990). 

For  the  presented  study,  the  VAX-Versions  of  “SF’  -  and  the  ITPP-algorithm  were 
used  (Klaes  et  al.,  1991;  Woolf,  1990).  “3F’  comparisons  were  made  with  LMD’s  IBM- 
Versions.  For  display  the  LMD-display  software  was  used  (Husson  et  al.,  1991). 

3.1  The  Retrieval  Algorithm 

The  “3F’  algorithm  (Chedin  et  al.,  1984;  Chedin  and  Scott,  1985;  Chedin  et  al., 
1985)  is  a  direct  physico-statistical  method.  For  the  retrieval  the  data  of  the  HlRS/2  and 
MSU  instruments  of  the  TOYS  instrument  package  are  used,  thus  representing  an  infiared 
and  a  mictowave  contribution. 

3.1.1  Theory 

The  philosophy  of  the  “3F’-method  is  threefold: 

•  The  method  must  be  physical,  thus  taking  into  account  all  information  of  the 
environment  and  the  observing  conditions  such  as  viewing  angle,  surface 
properties  like  emissivities,  land/sea  ^stinction,  surface  elevation,  solar  zenith 
angle,  etc.; 

•  The  algorithm  has  to  be  fast  because  of  the  large  amount  of  data  to  be  handled, 
dius  the  number  of  iterations  necessary  should  be  minimized; 
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The  first  guess  must  be  the  nearest  possible  to  the  actual  situation,  thus 
containing  as  much  a  priori  information  as  possible. 


♦ 


This  is  resolved  by  obtaining  the  first  guess  solution  from  the  observations 
themselves.  The  procedure  used  for  this  is  a  pattern  recognition  algorithm,  which  uses  die 
TOYS  Initial  Guess  Retrieval  (TIGR)  dataset  The  TIGR  dataset  is  composed  of  confuted 
sets  of  brighmess  temperatures  for  selected  atmospheric  conditions  at  40  levels  for 
(presently)  five  air  mass  types,  10  viewing  angles  and  different  ground  properties.  Since 
summer  1991  a  considerably  improved  TIGR-2  dataset  has  been  in  use  consisting  of  1761 
atmospheric  situations,  selected  out  of  more  tiian  150,000  radiosonde  observations. 

Transmissions,  radiances  and  weighting  functions  for  aU  of  the  channels  of  the 
instruments  used  in  the  “BP’-algorithm  and  for  tiie  different  satellites  (at  present  NOAA-8  to 
NOAA-ll)  are  calculated  by  u^g  die  Automated  Atmospheric  Absorption  Adas  (4A)  line- 
by-line  radiance  model  (Scott  and  Qiedin,  1981).  For  each  profile,  cloud-fiee  conditions 
are  assumed.  There  are  10  viewing  angles  between  0  and  60  deg.  (nadir  and  maximum 
scanning  angle  for  the  used  instruments,  respectively),  two  types  of  earth  surfaces  -  land 
and  sea,  and  ten  surface  pressure  values.  The  simulated  radiances  are  available  on  the  40 
“4A”  levels  (Table  4). 

In  the  proximity  recognition  algorithm  a  comparison  is  made  of  the  measured 
radiances  (by  the  satellite  instruments)  widi  die  archived  TIGR  radiances.  The  distance  d  in 
the  phase  space  of  the  TIGR-dataset  between  the  observed  brightness  temperatures  To  and 
the  brightness  temperatures  Tj  for  the  jth  TIGR  set  (j  =  1,...,1761),  is  given  for  a  selected 
number  of  TOVS-channels  n  as 


(3.1) 


where  Sk  is  the  variance  of  T/k)  over  “n”  TIGR  datasets.  The  TIGR  profile  with  the 
smallest  d  or  a  mean  of  a  number  of  profiles  around  die  most  proximate  situation  is  takoi  as 
die  initial  guess  for  die  rebievaL 

Since  the  TIGR-dataset  is  composed  of  cloud-free  conditions,  and  on  the  other 
hand,  many  of  the  HIRS/2  (infrared)  radiances  are  contaminated  by  clouds,  the  cloud 
influence  has  to  be  eliminated  before  the  retrieval  is  carried  out 
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To  detect  possible  cloud  contamination,  tests  are  carried  out  at  HIRS/2  sensor 
resolution.  These  tests  are  based  on  the  method  of  McMillin  et  al.  (1982),  which  has  been 
considerably  modified  (Wahiche  et  al.,  1986  and  Chedin,  1989). 

The  cloud  “clearing”  is  done  by  the  so-called  ^-method  (Chedin  et  al.,  1985  and 
Chedin,  1989).  The  ^-method  computes  a  clear  radiance  for  a  HIRS/2-channel  from  the 
difference  between  die  observed  and  the  TIGR-2  values  of  an  MSU  channel  which  has  the 
peak  of  its  energy  contribudon  near  the  one  of  die  HIRS/2-channel  under  consideration.  If 
the  microwave  radiances  are  not  contaminated  by  heavy  rainfall,  the  integration  of  the 
energy  contribution  over  approximately  the  same  part  of  the  atmosphere  yields  a  difference, 
which  describes  the  major  part  of  the  initial  value  and  the  equivalent  clear  value  of  the 
HIRS/2  channel.  See  Appendix  A  for  details  of  the  ^-method. 

Having  obtained  the  first  guess  soluticMi  the  temperature  retrieval  is  carried  out  The 
method  of  temperature  retrieval  is  based  on  a  Bayesian  approach  (Chedin  et  al.,  1985).  In 
principle  the  basis  is  the  equation  of  radiative  transfer  in  the  form 

TB(y,)  =  /  4>*{v„()r(()d(  (3.2) 

where  TB(v)  denote  the  satellite  measured  brightness  temperatures,  T  die  atmospheric  tem¬ 
peratures  to  be  retrieved,  is  the  kernel,  which  includes  the  atmo^heric  transmittances, 
^  is  the  altitude  variable  (arbitrary,  generally  the  atmospheric  pressure),  and  the  V{  denote 
the  various  observation  conditions  (instrument  channels,  viewing  angle,  surface 
characteristics,  etc.). 

Linearizing  Equation  3.2  by  quadrature  we  arrive  at 

TB  =  ^*T  +  €  (3.3) 


where  e  indicates  non  systematic  error.  For  the  errors  several  assumptions  are  made:  they 
should  have  a  mean  of  zero,  should  be  normal  distributed  in  terms  of  their  probability  and 
their  error  covariance  matnx  Se  should  be  known;  the  temperatures  T  are  assumed  to  be 
described  by  a  multivariate  normal  distribution  with  mean  Tq  and  covariance  mattix  5r< 
Applying  Bayes’  conditional  probahUity  rule 


'  '  P{TB) 


(3.4) 
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we  find  a  probability  density  function  of  the  atmospheric  state  T  under  the  assumption  that 
the  observations  are  given.  Finding  the  profile  (atmospheric  state)  with  the  maximum 
probability  we  obtain  an  estimator  -  nuudmum  a  posterior  (MAP)  -  as 

where  5T  =  -  To ;  STB  =  TB^  -  TBq,  with  Tq  and  TBq  the  expectation  values 

respectively.  The  initial  guess  retrieval  procedure  has  no  bias  (Chedin  et  al.,  1985),  thus 
the  means  ofT*  =  Tuap  -  Tfo  and  TB*  =  TBg  -  TBjg  are  zero.  [X]'l  denotes  the  inverse 
and  the  transpose  of  a  matrix  X.  It  follows 

5r*=  (3-6) 

dTB 

The  matrix  4>*  contains  the  Jacoblans - i  and  i  are  the  levels. 

dTj 

3.1.2  Application 

The  retrieval  is  carried  out  on  so-called  “boxes”  of  1(X)  km  x  1(X)  km,  containing 
four  to  six  HIRS-2  spots  which  are  merged  into  the  MSU  spots  (see  Figure  4),  taking  into 
account  the  different  ground  resolution  of  the  TOYS  instruments  used  for  the  retrieval. 

After  having  completed  the  temperature  retrieval,  cloud  parameters  are  estimated  by 
the  “method  of  coherence  affective  cloud  amounts”  (Wahiche  et  al.,  1985),  a  con^lex 
procedure  yielding  cloud  top  pressure,  cloud  top  temperature  and  effective  cloud  amount 
(cloud  cover  times  cloud  emissivity). 

The  extraction  of  water  vapor  information  out  of  the  TOYS  observations  is  done 
after  the  temperature  profiles  have  been  estimated.  The  inversion  of  water  vapor  para¬ 
meters  is  done  for  the  three  atmospheric  layers  300-500  hPa,  500-800  hPa  and  800- 
1000  hPa.  The  parameters  retrieved  include  water  vapor  contents  for  the  tiuee  layers,  total 
prec^itable  watCT  and  surface  tenq)eraturB. 

The  method  used  ha«  is  a  ridge-type  approximation  for  a  simultaneous  retrieval  of 
the  above  mentioned  parameters.  Assuming  that  Sj  =  Ojl  for  the  HIRS/2  channels  8, 10, 
11  and  12,  and  further  5^  =  a]!,  where  I  is  the  unity  matrix  we  arrive  at  a  ridge-type 
estimator  (Chedin  et  al.,  1985) 

5^;  =  [  0  *^0  *+//]■'[  0*]  STB  *  (3.7) 
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for  q*  =  Ui,i  =  for  the  three  layers,  f/,  is  the  relative  humidity,  =  7^.  The  8 
denote  the  deviation  firom  the  fost  guess. 

After  all  these  operations  have  been  carried  out  further  derived  geophysical 
parameters  can  be  estimated,  e.g.  geopotential  thicknesses  or  thermal  winds,  by  using  the 
retrieved  temperatures  and  humidity  parameters  to  derive  the  virtual  temperatures  for 
chosen  layers  and  apply  the  it.>ults  to  the  reflective  equations. 

The  “ST’  algorithm  can  be  inproved  by  providing  additional  information  in  the  fonn 
of  forecast  fields  for  the  refinement  of  the  initial  guess.  The  final  retrieval  is  carried  out 
with  die  thus  specified  initial  guess. 

3.2  The  ITPP  Retrieval  Algorithm 

The  “International  TOYS  Processing  Package”  (TTPP),  also  known  as  the  ‘TOYS 
Export  Package”,  has  been  developed  at  the  Cooperative  Institute  for  Meteorological 
Satellite  Studies  (OMSS)  at  the  University  of  Wisconsin,  Madison,  Wisconsin  (Smith  et 
al.,  1985).  The  first  physical  approach  consisted  of  the  iterative  scheme  developed  by 
Smith  (1968).  The  current  version  permits  the  simultaneous  retrieval  of  surface 
temperature,  temperature  and  moisture  profiles  (Smith  et  al.,  1985).  The  physical  character 
makes  it  possible  to  consider  the  influence  of  surface  parameters  (elevation,  emissivity, 
etc.)  and  of  clouds. 

The  essential  retrieval  method  is  the  simultaneous  temperature  and  water  vapor 
retrieval  (Smith  and  Woolf,  1984).  The  basis  again  is  as  in  (3.2),  the  equation  of  radiative 
transfer  in  the  form 

ra,  =j4>{v„p)r(p)<4>  (3.8) 

where  the  pressure  p  rq)laced  the  vertical  coordinate  ^  of  (3.2),  TBi  are  the  radiance 
expressed  as  brightness  temperatures  as  observed  by  the  satellite,  ^Vip)  is  the  kernel  and 
contains  a  pressure  dependency  of  the  profiles  to  retrieve,  T(p)  is  the  atmospheric 
temperature  profile  to  be  retrieved. 

Integrating  the  radiative  transfer  equation  by  parts  and  developing  a  small 
poturbation  of  the  brightness  temperature  we  obtain 
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where  &c  =  x-xq,  and  xo  is  an  a  priori  mean  of  first  guess  of  the  variable  x. 

Each  perturbation  is  assumed  as  dependent  on  a  pressure  function  ^p)  thus  that 

&l{p)  =  ^MPpiip)  (3. 10) 

i=l 


and 


^(p)  =  - 


i=W+l 


(3.11) 


with  q(p)  being  the  water  vapor  mixing  ratio  and  g  the  gravity  acceleration.  With  the 
hydrostatic  equation  and  the  gas  law  Equation  3.10  becomes 


i=i 


(3.12) 


Setting  ao  =  ST s  and  substituting  Equations  3.11  and  3.12  into  3.9  yields  the 
following  matrix  equation 

=  (3.13) 

1*0 


where  K  denotes  the  number  of  sounding  channels  used  for  the  retrieval.  Thus, 

STB]  =  ditt.  (3.14) 

If  the  number  of  channels  K  is  greater  or  equal  to  the  number  of  variables  plus  one 
(M+1 ),  then  least  square  methods  transform  3. 14  into 


a  =  * 


(3.15) 


which  can  be  approximately  written  as 


a  =  +  //]■  ♦  (3.16) 

The  term  with  the  Lagrange  parameter  /(nominally  0.1)  (Smith  et  al.,  1985)  is 
introduced  to  stabilize  the  matrix  inversion.  Again,  /  denotes  the  unity  matrix,  [A]-^  the 
inverse  and  [Af  the  transpose  of  matrix  A. 

The  choice  of  the  pressure  functions  ^(p)  is  arbitrary;  however,  in  the  ITPP  the 
profile  weighting  functions  of  the  equation  of  radiative  transfer  are  used.  One  other 
possibility  was  the  use  of  eigenvectors  of  the  temperature  and  water  vapor  covariance 
matrices. 

The  initial  profile  in  the  processing  of  the  TOYS  retrieval  may  be  obtained  from 
climatology  or  a  regression  based  on  synthetic  radiances  or  on  radiances  obtained  by 
matching  radiosonde  observations.  Both  possibilities  are  part  of  the  ITPP.  A  third  pos¬ 
sibility  is  the  introduction  of  the  first  guess  profiles  fit>m  analysis  or  forecast  information  - 
this  information  has  to  be  provided  by  the  user  -  from  outside. 

The  HIRS/2  and  MSU  radiances  are  used  on  a  3  x  3  field  of  HIRS/2  field-of-views 
(FOV),  which  is  about  75  km  resolution  (Smith  et  al.,  1985).  The  cloud  contamination  is 
treated  by  the  iV*  method. 

3.3  Pre-Processing 

For  both  retrieval  packages  the  ITPP  pre-processing  (Woolf,  1990)  was  applied. 
This  includes  1)  the  treatment  of  the  direct  ingested  raw  data  of  the  HRPT  (High 
Resolution  Picture  Transmission)  data  stream,  2)  the  extraction  of  the  TIP  (TIROS 
Information  Processor)  data,  which  include  the  TOYS  data,  3)  the  extraction  of  the  TOYS 
data,  4)  the  calibration  and  earth-location  of  the  HIRS/2  and  MSU  data  (see  e.g.  Planet, 
1988  for  more  details  on  the  data  format),  and  5)  the  limb-correction,  depending  on 
whether  “37”  or  ITPP  is  used.  Thus,  both  packages  use  the  same  raw  data.  (This  is  also 
an  advantage  for  the  operational  tqrplication  of  both  packages.) 
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Chapter  4 


Weather  Situation 

For  the  application  of  the  two  retrieval  schemes  the  weather  situation  for  13  August 
1990  over  Europe  was  chosen.  On  one  hand  there  has  been  the  mesoscale  phenomenon 
which  has  to  be  investigated,  a  cloud  cluster  over  eastern  France  and  western  Germany  in 
the  synoptic  surroundings  of  a  decreasing  high  pressure  situation.  On  the  other  hand  there 
has  been  the  direct  ingested  HRFT  data  available,  received  directly  at  the  German  Military 
Geophysical  Office  (GMGO)  at  Traben-Trarbach,  Germany  (Klaes,  1990).  The  TIP-data 
were  extracted  out  of  the  HRPT-data  and  preprocessed  to  be  ready  as  “37”  as  well  as  ITPP 
input  data. 

The  synoptic  situation  for  13  August  1990  over  Europe  is  characterized  by  hot 
subtropical  air  over  Southern  and  Central  Europe.  Frequently  the  maximum  temperatures 
over  the  continent  reached  more  than  30  deg  C.  The  synoptic  surface  pressure  field  was 
characterized  by  high  pressure  over  the  Karpates  and  southern  France/Italy  (Figure  5). 
Combined  with  this  surface  field  there  was  an  upper  ridge  over  Central  Europe  which 
extended  at  00  UTC  from  die  western  Baltic  Sea  to  Southern  France  (Hgure  6). 

Two  weak  upper  troughs  were  analyzed  over  The  Netherlands  and  over  the  Gulf  of 
Biscay;  the  latter  being  the  stronger  one.  The  eastward  movement  of  the  upper  ridge  as 
well  as  the  strong  heating  during  the  day  caured  a  wealrening  of  the  high  pressure  influeice 
over  Europe.  There  was,  according  to  the  upper  trough  and  the  heating  of  the  surface,  a 
flat  low  pressure  system  over  The  Netherlands,  with  a  frontal  system  along  the  North  Sea 
coast  and  over  the  Gulf  of  Biscay,  mainly  recognizable  at  upper  levels. 

A  mature  low  pressure  system  is  nearly  stationary  west  of  Iceland,  its  already 
occluded  frontal  system  extending  from  the  North  Atlantic  over  Northern  Scandinavia,  the 
North  Sea  and  the  British  Isles.  The  fronts  ate  already  weak. 

Thus,  the  mesoscale  situation  is  not  the  ideal  case  -  ideal  case  meaning  that  the 
mesoscale  phenomenon  is  situated  in  a  homogeneous  synoptic  environment  -  a  situation 
which  will  be  found  very  rarely  over  Central  Europe.  The  present  situation  is  characterized 
by  unstable,  subtropical  air  over  Central  Europe,  being  further  destabilized  by  strong 
heating  as  well  as  an  approaching  upper-air  trough. 

The  satellite  picture  (Hgure  7)  is  an  AVHRR  Channel  4  (10.S  pm)  composite  of  the 
three  orbits  9704, 9705  and  9706  of  NOAA-11  over  central  Europe,  projected  on  a  polar 
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stereographic  projection  by  the  Free  University  of  Berlin  (1990).  The  structures  discussed 
above  can  be  easily  detected.  Since  the  time  of  the  displayed  orbits  is  about  12  hours  later 
than  the  synoptic  maps  shown,  there  are  a  lot  of  thunderstorm  cells  developed  along  the 
frontal  zone.  Two  structures  are  most  interesting:  the  big  cloud  cluster  over  eastern 
France/south  western  Germany,  which  has  a  horizontal  extension  of  about  300  x  100  km^ 
and  a  circular  structure  over  the  Gulf  of  Biscay,  also  a  well  developed  structure  of  about 
100  X  100  km2 

The  two  retrieval  algorithms  discussed  in  Chapter  2  were  applied  to  direct  readout 
satellite  data  for  all  the  cases  discussed  above.  In  the  following  chapter  some  of  the  results 
will  be  shown  and  some  aspects  of  the  results  will  be  discussed  in  the  light  of  the 
qualitative  comparison  with  “conventional'*  meteorological  data  and  analyses  discussed 
above. 
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Chapter  5 


Retrieval  Results 

In  the  following,  results  of  the  two  retrieval  procedures  “31-2”  and  ITPP  V4.1  will 
be  shown.  Both  versions  were  available  and  used  as  VAX  versions  (Klaes  et  al.,  1991a, 
Woolf,  1990).  “3r’  was  also  used  in  the  original  IBM  version  at  the  C.I.R.C.E.  (Centre 
Inter  Regional  de  Calcul  Electronique  du  C.N.R.S.)  at  Orsay,  France. 

The  VAX  versions  were  installed  on  the  AIMS  (Air  Force  Interactive 
Meteorological  System)  computer  system  of  the  Geophysics  Directorate,  Satellite  Meteor¬ 
ology  Branch,  Hanscom  Air  Force  Base,  Massachusetts,  U.S.A.  (Klaes,  1991). 
Comparisons  of  the  “3r’  VAX  version  have  been  done  with  the  original  version  of  “3r’  at 
Palaiseau.  The  results  were  found  identical  (Klaes  et  aL,  1991a). 

Comparisons  have  been  done  between  “31-2”  and  “31-1”  for  the  case  investigated 
here,  the  results  have  been  discussed  elsewhere  (Klaes  et  al.,  1991b). 

The  results  displayed  in  this  chapter  are  geophysical  fields  of  both  retrieval 
packages.  There  are  examples  of  the  Geopotraitial  Thicknesses  500/1000  hPa  as  a  measure 
of  the  average  temperature  of  the  lower  half  of  the  atmosphere  as  well  as  the  quantity  total 
precipitable  water  as  a  measure  of  the  water  vtqwr  contents  of  the  atmosphere.  There  are 
some  conq>arisons  of  vertical  profiles  of  tonperature,  they  are  compared  with  radiosonde 
measurements  of  selected  places.  The  data  used  are  the  NOAA-11  orbits  of  13  August 
1990,  1233  UTC  and  1403  UTC.  The  geophysical  fields  retrieved  from  them  are 
displayed  together  and  therefore  are  referred  to  as  1315  UTC. 

5.1  “31”  Results 

Figure  8  shows  the  Qear/Not  Clear  Flags  of  the  two  combined  orbits  of  1315 
UTC,  13  August  1990.  The  results  are  displayed  for  the  “3?’  “boxes”  (cf.  Figure  4).  The 
green  spots  denote  cloud-free  boxes  over  land,  the  blue  spots  denote  cloud-free  boxes  over 
sea,  the  yellowish  boxes  denote  cloud  contaminated  boxes.  The  two  cloud  cluster  areas  are 
clearly  represented,  the  one  over  the  Gulf  of  Biscay  more  clearly  because  of  its  isolated 
nature.  The  one  over  central  Europe  is  embedded  in  the  other  clouds  of  the  frontal  zone. 
Note  that  even  the  area  of  low  cumulus  cloi^  over  Eastern  Europe  is  recognized  as  cloudy 
(conq>are  also  the  AVHRR  image  Figure  7).  Looking  at  the  cloud  top  pressures  obtained 
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after  the  temperature  retrieval  (Figure  9)  it  can  be  seen  that  tte  regions  of  the  cloud  clusters 
are  detected  as  regions  with  high  reaching  clouds  -  cloud  top  pressures  of  350-200  hPa. 
Nevertheless,  this  is  indicated  by  black  spots,  a  number  of  boxes  is  rejected  because  of  too 
large  cloud  amounts,  thus  no  estimation  of  cloud  altitude  possible.  Also  apparent  is  that  the 
frontal  area  over  the  British  Isles  is  composnl  of  lower  clouds. 

Figure  10  shows  the  geopotential  thicknesses  S(X)/1(XX)  hPa  as  point  plots  for  the 
1315  UTC  date.  It  can  be  seen  that  the  frontal  zone  over  the  North  Sea/British  Isles  is 
represented  as  a  strong  gradient  in  the  thickness  field.  The  mature  low  west  of  Iceland  is 
represented  as  a  center  of  cold  air  Gow  thicknesses). 

The  area  of  the  mesoscale  cloud  cluster  over  Eastern  Europe,  southwestern 
Germany  is  represented  as  a  cyclonic  deformation  (trough)  in  the  thickness  field  exactly 
over  the  area  (compare  the  AVHRR  image  Figure  7).  The  warm  subtropical  air  approach¬ 
ing  western  Europe  from  north  Africa  can  be  seen  as  a  tongue  of  warm  air  from  Maroc. 
The  results  are  confirmed  if  the  field  is  displayed  as  an  isoline  field  (Figure  1 1).  Results 
have  to  be  looked  at  first  in  terms  of  boxes  before  they  are  displayed  as  isoline  fields  to 
preclude  spurious  contouring  errors. 

The  moisture  field  in  terms  of  total  wato:  vapor  content  (total  precipitable  water)  is 
displayed  in  Figures  12  and  13.  There  is  a  well  expressed  moist  maximum  in  the  zone  of 
the  frontal  area  over  Scandinavia  and  the  British  Isles,  the  moisture  field  confirms  that  the 
front  is  weak  in  its  middle  part  over  the  North  Sea.  The  dry  zone  over  Central  Europe 
denotes  very  significantly  the  area  of  dry  air  which  brought  the  longest  period  of  drought 
over  central  Europe  since  the  dry  summer  of  1976  (Deutscher  Wetterdienst,  1990).  It  can 
also  be  seen  that  the  subtropical  air  approaching  from  the  southwest  is  carrying  more 
moisture  witii  it 

5.2  “ITPP”  Results 

Displaying  the  geopotential  thicknesses  500/1000  hPa  of  13.15  UTC  of  the 
13  August  1991  as  point  plots  of  the  ITPP  V4.1  (Figure  14)  we  see  that  the  stractures  of 
the  temperature  field  are  retrieved  in  a  similar  way  -  the  cloud  cluster  area  over  eastern 
France  and  southwestern  Germany  is  also  re|n:esented  as  a  deformation  in  the  thickness 
gradient  which  may  be  easily  seen  in  the  isoline  display  (Figure  15). 

The  moisture  field  in  terms  of  total  precipitable  water  (Figures  16  and  17,  point 
plots  and  isoline  displays  respectively)  shows  the  frontal  structure  as  an  increase  of 


moisture.  The  dry  area  over  central  Europe  is  smaller.  Generally  the  amounts  of  ITPP 
retrieved  moisture  are  higher  than  with  “ST*. 

Statistical  results  of  the  comparison  as  well  as  comparisons  with  radiosondes  are 
contained  in  Table  5. 

5.3  Vertical  Profiles 

Hgure  18  shows  an  example  of  a  vertical  profile  of  temperature  obtained  by  the  two 
retrieval  algorithms  “ST*  and  ITPP  as  well  as  by  a  radiosonde  observation.  The  v^cal 
sounding  of  box  8  of  HIRS  line  68  is  chosen,  being  situated  nearest  to  the  radiosonde 
10338  Hannover  in  Germany.  Displayed  are  the  results  of  “31-2”  (solid  thick),  ITPP  V4 
(dashed)  and  the  radiosonde  (thin  solid).  It  can  be  seen  that  in  the  region  between  900  hPa 
and  300  hPa  the  vertical  stability  of  the  atmosphere  is  simulated  similar  to  the  radiosonde 
from  both  retrieval  algorithms.  This  can  be  ^own  for  all  the  retrievals  near  radiosonde 
stations.  It  can,  therefore,  be  concluded  that  the  thermal  stability  of  die  troposphere  is  weU 
retrieved  by  eidier  algorithm. 

There  is  a  difficulty  in  the  treatment  of  the  surface.  This  is  explainable  as:  1)  one  of 
the  major  difficulties  in  retrieval;  2)  the  fact  that  the  day  under  investigation  was  a  very  hot 
day  in  terms  of  surface  temperature  -  in  most  parts  of  Gormany  the  temperatures  reached  or 
exceeded  303  K;  3)  that  there  is  a  difference  in  the  terrain  elevation  of  the  radiosonde 
station  and  the  mean  terrain  altitude  of  the  box  or  die  retrieval  area  respectively  (composed 
of  4-6  HIRS/2  spots).  Thus,  there  may  be  a  major  difference  in  surface  temperatures.  If 
the  tenqietatures  at  the  same  levels  are  compared,  we  obtain  a  difference  of  approximately  3 
K  at  the  surface. 

Note  that  the  character  of  the  retrieved  profiles  and  radiosoundings  are  essentially 
different  The  tadiosoundings  are  taken  over  a  time  of  about  100  minutes  for  an  area  op  to 
400  km  in  diameter,  depending  on  the  wind  drift  of  the  sonde,  whereas  retrievals  are  the 
results  of  instantaneous  measurements,  representative  for  one  moment  the  moment  of  the 
radiance  measurement  of  the  satellite  instrument  valid  for  the  area  of  the  resolution  of  the 
ground  spot  -  in  case  of  more  than  one  instrument  the  area  of  the  coarsest  resolution 
instrument 

The  radiosonde  is  making  a  point  measurement  at  the  position  in  space  where  the 
sonde  is  during  the  ascent  Retrieved  profiles  represent  volume  avoages,  the  volume  being 
defined  horizontally  by  the  size  of  the  footyrint  of  the  respective  instrument  and  in  the 
vertical  by  the  widtii  of  the  respective  contribution  (weighting)  function. 
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Thus,  radiosonde  measurements  and  retrievals  must  not  be  compared  without 
taking  into  consideration  the  different  character  of  both  measurements.  The  comparison  of 
radiosonde  measurements  with  retrieval  results  of  13  August  1990  is  based  on  four 
overpasses  of  NOAA-10  and  NOAA-11  spacecraft  The  collocation  criteria  were  the 
following:  1)  the  collocated  retrieval  profiles  must  be  located  within  100  km  diameter 
around  the  radiosonde  observation;  and  2)  the  time  difference  must  be  less  than  3  hours. 
Given  these  conditions,  534  collocated  retrievals  could  be  found  to  compare  the  thicknesses 
500/1000.  Based  on  the  mean  of  5535.17  gpm  the  bias  was  found  to  be  -34.65  gpm  (“3F’ 
retrieval  minus  radiosonde  measurement),  corresponding  to  -1.71  K.  The  standard 
deviation  is  20.39  gpm  (corresponding  to  1.01  K). 

Total  water  content  was  found  to  be  biased  by  -0.22  cm  with  a  standard  deviation 
of  0.57  cm,  based  on  the  radiosonde  average  of  2.16  cm.  The  comparison  results  are 
summarized  in  Table  5. 

The  temperature  measurements  in  terms  of  thicknesses  can  be  considered  to  be  of 
comparable  range  to  each  other,  the  retrieval  bemg  colder  than  the  collocated  radiosonde 
profiles.  The  humidity  in  terms  of  total  water  vapor  contents  is  underestimated  by  the 
satellite  retrieval,  with  an  error  of  about  26%.  The  majority  of  this  di^erence  is  due  to  the 
underestimation  of  the  humidity  in  the  lowest  retrieval  level,  where  the  highest  amount  of 
humidity  occurs,  which  cannot  be  resolved  by  the  HERS  instmment  channel  10. 

It  can  be  stated  in  general  that  the  retrieved  profiles  are  in  accqitable  agreement  with 
the  in  situ  measurements. 
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Chapter  6 


Conclusions 

The  purpose  of  this  reserach  is  to  demonstrate  the  ability  and  capability  of  two 
TOYS  retrieval  algorithms,  “Sr*  and  the  ITPP. 

It  is  demonstrated  how  the  two  retrieval  algorithms  treat  the  temperature  retrieval  of 
TOYS  instrument  soundings  as  well  as  the  deduction  of  water  vapor  quantities.  It  is 
shown  that  the  retrieval  algorithms  are  able  to  produce  geophysical  fields  over  a  specified 
area,  according  to  the  satellite  passes.  Meteorological  features  like  mesoscale  cloud 
clusters,  down  to  a  resolution  of  the  meso-^  scale  are  resolved.  It  must  be  emphasized  that 
for  this  study  TOYS  data  alone  were  used  for  the  retrievals.  No  forecasts  were  used  as 
additional  information. 

It  is  demonstrated  from  Figures  19  and  20  that  the  derived  fields  fit  well  the  results 
obtained  by  “conventional”  analysis.  Figure  19  shows  the  geopotential  thicknesses 
500/1000  hPa,  analyzed  by  the  German  Weather  Service  for  OOUTC,  13  August  1990. 
The  analysis  is  based  on  all  available  observations,  including  satellite  data.  Figure  20 
shows  again  the  result  of  “BF’  for  1315  UTC  of  the  same  day.  Considering  that  there  is 
more  than  12  hours  difference  between  the  two  fields,  it  can  be  concluded  that  there  is 
considerable  agreement  between  the  two.  Again,  it  is  emphasized  that  the  mesoscale  cloud 
cluster  as  a  deformation  in  the  thickness  field  is  recognized  properly  by  the  retrieval. 

Additionally  it  can  be  stated  that  the  cold  air  regions  in  the  area  of  the  mature  low 
pressure  system  west  of  Iceland  as  well  as  those  over  Sweden/Finland  and  western  Russia 
are  well  simulated.  It  is  concluded  that  the  retrieval  algorithms  are  a  good  instrument  to 
derive  geophysical  fields  from  satellite  data  alone,  especially  in  between  synoptic 
observation  times  and/or  in  regions  whoe  “conventional”  observations  are  sparse.  This  is 
supported  by  the  good  simulation  of  the  vertical  stabili^  of  the  troposphere  in  the  layer  900 
to  BOOhPa. 

The  humidity  retrievals  are  not  that  well  defined;  it  is  not  possible  from  tiiis  study 
alone  to  state  accuracies  (as  well  we  do  not  know  the  “truth”).  Comparisons  need  to  be 
made  with  radiosonde  data  as  well  as  with  data  from  different  satellite  sensors  such  as  the 
SSM/I  and  SSM/T-2  (microwave)  derived  moisture  fields. 

To  make  statistical  statements  additional  research  is  required  with  more  data  - 
especially  radiosonde  and  passive  microwave  satellite  data. 
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Hgure  4:  “SF’  BOXES:  HIRS/2  spots  merged  into  MSU  spots. 
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Figure  5:  Surface  pressure,  analyzed  by  Deutscher  Wetterdienst  for  13  August  1990, 
OOUTC 


Figure  6:  Geopotential  height  500  hPa,  analyzed  by  Deutscher  Wetterdienst  for 
13  August  1990,  00  UTC. 
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Figure  7:  AVHRR  channel  4  (10.5^iin)coniposite  image  of  NOAA-11  orbits  9704, 
9705  and  9706  for  13  August  1990  (central  time  13.15  UTC),  processed  by 
the  Free  University  of  Berlin  (1990). 


Hgore  11:  Geq)Otential  thicknesses  500/1000  hPa  /gpdam,  retrieved  by  “31-2”  for 
13  August  1990  at  13.15  UTC  from  NOAA-11  orbits  9705  and  9706. 

Results  as  before  but  as  isolines. 
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Figure  12:  Total  water  vapor  contents  as  PPW  TOT  /cm,  retrieved  by  “31-2”  for 


13  August  1190  at  13.15  UTC  from  NOAA-11  orbits  9705  and  9706. 
Results  of  “boxes”. 


Figure  13:  Total  water  contents  as  PPW  TOT  /cm,  retrieved  by  “31-2”  for  13  August 


1990  at  13.15  UTC  from  NOAA-11  orbits  9705  and  9706.  Results  as  for 
Figure  12  but  as  isolines. 
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Hgure  14:  Geopotential  thicknesses  500/1000  hPa  /gpdam,  retrieved  by  “TTPP  4”  for 
13  August  1990  at  13.15  UTC  from  NOAA-11  oifrits  9705  and  9706. 
Results  of  75  km  areas. 


Figure  15;  Geopotential  thicknesses  500/1000  hPa  /gpdam,  retrieved  by  “ITPP  4”  for 
13  August  1990  at  13.15  UTC  from  NOAA-11  orbits  9705  and  9706. 


Results  as  for  Hgure  14  but  as  isolines. 
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SAT:  NOAA  11 
DATE . :  130890 


UVT...  52  93  TIME  (UT)  :  121818 

LON..  9  80  CLEAR -SPOTS/%.  00 


Figure  18:  Vertical  profiles  of  atmospheric  temperature,  obtained  by  “31”  (thick  solid 
line),  ITPP  (dashed  line)  for  13  August  1990  at  13.15  UTC  from  NOAA-1 1 
orbits  9705  and  9706  and  obtained  from  the  radiosonde  observation  at 
Hannover  (10338)  on  13  August  1990, 1200  UTC. 
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Figure  19;  Gcopotential  thicknesses  500/1000  hPa,  analyzed  by  Deutscher  Wetterdienst  for  13  August  1990,  00  UTC. 
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Rgure  20:  Geopotential  thicknesses  500/1000  hPa  /gpdam,  retrieved  by  “31-2”  for 
13  August  1990  at  13.15  UTC  from  NOAA-11  orbits  9705  and  9706. 
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Table  1:  Qiaiacteristics  of  the  HIRS/2  sounding  channels. 
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Table  2:  Characteristics  of  the  MSU  sounding  channels. 


MSU  CHARACTERISTICS 

MSU 

Channel 

^^Numb^_ 

1 

55.73 

^^^dow 

Surface 

2 

53.73 

700 

3 

54.96 

02 

300 

4 

57.95 

Oi 

90 

Table  3:  Characteristics  of  die  SSU  sounding  channels. 


SSU  CHARACTERISTICS 

SSU 

Channel 

Number 

Wavelength 

/Mni 

Principal 

Absorbing 

Ccmstituents 

Level  of 

Peak  Energy 
CcHitribution/iiPa 

1 

WSM 

CQz 

15.00 

2 

002 

4.00 

3 

15.00 

002 

1.50 
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Table  S:  Statistics  (4  oibits)  “31  Std”  OMnpared  with  values  derived  from  radiosondes, 
13  August  1990. 

Mean 


Thicknesses 

Radiosondes 
(gpm)  _ 

Bias  (gpm) 
or(K) 

STOV  (gpm) 
or  (K) 

#  Items 

socviooo 

5535.17 

-34.65 

20.39 

534 

-1.71 

1.01 

Total  water 

(cm) 

(cm) 

(cm) 

#itms 

4iot 

2.16 

-0.22 

0.57 

437 
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Appendix  A 


The  ^-Method  in  “31” 


The  ^-Method  uses  the  non  sensitivity  of  the  MSU  channels  to  clouds  (provided  there  is 
not  too  much  rain  coitaminadcm)  and  tries  to  create  “cleared”  HIRS/2  channel  radiances  by 
comparing  MSU  brightness  temperatures  with  HIRS  brightness  temperatures  for  charuiels, 
whose  energy  omtribution  functions  peaks  at  the  same  level 
Thus  Y-radiances  are  created  as 

*''“(<•) (A.1) 

where  i  =  4, 5, 6, 15  while  j=  3, 2, 2, 2,  respectively,  and 

•P““(14)  =  ho  +  S  h.T""®(0  +  ^  CjT^^^'U)  (A.2) 

‘  J 

where  i  =  and  15;  for  j  =  2,3. 
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